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Consider a discrete-time quantum walk on the iV-cycle governed by the following condition: at 



J2 
0,: 

every time step of the walk, the option persists, with probability p, of exercising a projective 
^^| measurement on the coin degree of freedom. For a bipartite quantum system of this kind, 

we prove that the von Neumann entropy of the total density operator converges to its maxi- 
mum value. Thus, when influenced by decoherence, the mutual information between the two 
subsystems, corresponding respectively to the space of the coin and the space of the walker, 
eventually must diminish to zero. To put it plainly, any level of decoherence greater than zero 
VJ ■ forces the system eventually to become completely "disentangled" . 
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1. Introduction 

A quantum walk (QW) is a reversible process commonly described as the quantum- 
mechanical analogue of a classical random walk. In recent times, quantum walks have at- 
tracted extensive attention, mainly for their value as potential sources of new algorithms 
flKempe 2003] |Ajnbainis^003| iKendon 2006| [Venegas-Andraca 2008| [ Konno 2008J) . 
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Like classical random walks, quantum walks are classified into two main types: discrete- 
time quantum walks QNayak and Vishwanath 2000 lAmbainis et al. 2001IIAharanov et al. 2001|) 



and continuous-time quantum walks (jFahri and Gutmann 19 98; Chil ds et al. 2002)) . Both 
types exhibit similar dynamical properties, but the continuous type can be obtained from 
the discrete type by way of a suitable limiting process (|Strauch 2006)IChilds 2010]) . 

Whereas continuous-time quantum walks can be modeled in terms of a single position 
Hilbert space Tip, discrete-time quantum walks cannot be modeled except in terms of 
a bipartite system involving the tensor product of the position Hilbert space Tip and 
an auxiliary coin Hilbert space T-Lq- The conditional shift operator, which governs the 
itinerary of the walker, induces entanglement between the degree of freedom of the coin 
and the spatial degree of freedom of the walker. 

To capture the full computational power inherent in any quantum-mechanical pro- 
cess, it is essential to be able to model the phenomenon of quantum entanglement. Re- 
cently, several important studies of quantum entanglement in discrete-time quantum 
walks have been completed, both from a numerical perspective (jCarneiro et al. 2005| 
Mal oyer and Kendon 2007||Venegas-Andraca and Bose 2009[ ) and from an analytical per- 
spective (|Abal et al. 2006) lAnnabestani et al. 2010a|) . 

In recent years, several schemes have been proposed to implement quantum walks 
in realistic media ( |Travaglione and Milburn 2002| IDur et al. 20021 ISanders et al. 20031 
IDu et al. 2 0031 |Ryan et al. 2005[lEckert et al. 2 00JllZcjieLaL_2Q061 [va"n Hoogdalem and Blaauboer 2009[ ). 
However, any attempt to implement a quantum system in a physical channel must take 
into consideration the critical issue of "decoherence" , whereby the idiosyncratic features 
of a quantum system succumb to macroscopic (classical) manifestations. Various mathe- 
matical models of decoherence in discrete-time QWs have been investigated both numeri- 
cally and analytically ( |Kendon and Tregenna 200~2 IBrun et al. 2003allBrun et al. 2003b[ 



|Kendon and Tregenna " 2003a; K endon and Tregenna 2 003b; R omanelli 2 005 Ko sik et al. 20061 



IRichter 20 07; Zhang 2008; Banerjee et al. 2008|lLiu~ahd Petulante 2010; Annabcsta ni et al. 2010b)) . 



In this paper, as in (Brun et al. 2003a; B run et al. 2003b)) . we adopt a specific model of 
decoherence, which might best be described as follows. At each time step of the quantum 
walk, an observer stands ready to perform a projective measurement on the coin degree 
of freedom. The probability of performing a measurement is given by a fixed parameter 
p, called the "decoherence rate" . 

In the current literature, not much evidence can be found of attempts to provide a 
precise formulation of the relationship between entanglement and decoherence. A notable 
exception is the numerical study given by ( |Maloyer and Kendon 2 007). However, the 
measure of entanglement employed in ( |Maloyer and Kendon 2007| ), called the negativity 
measure, differs considerably from the measure of entanglement adopted in this paper. 

Here, we utilize the concept of von Neumann entropy to quantify the information 
content of the various components of the quantum walk system, including the mutual 
information between its subsystems (coin and position). This approach enables us to 
provide a precise formulation of the measure of entanglement between the subsystems. 
In the presence of any non-zero level of decoherence, the von Neumann entropy associated 
with the total probability density of the system tends to its maximum value, implying 
the total collapse of entanglement between subsystems. 
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2. Quantum walks on the TV-cycle subject to decohrence 

We proceed to define the elements, as formulated in ()Liu and Pctulantc 2 010|) . of a quan- 
tum random walk on the TV-cycle. The definitions and corresponding notations are anal- 
ogous to those outlined in (jLiu and Petulante 20f0|) . Let t denote the number of time 
steps from the moment the discrete quantum walk is launched on the iV-cycle. The tem- 
poral evolution of the system is modeled by ipt — U tr ipo, where tpo is the initial state 
and U = S(I (?) Uq) is the evolution operator on the Hilbert space % = Hn ® %• Here 
Uq denotes the coin operator on the coin subspace %2 spanned by an orthonormal basis 
{\j)ij — 1>2.}, and S denotes the conditional shift operator on the position subspace 
T-Ln spanned by an orthonormal basis {|x),<E € Zjv}. Thus, a typical state ip in % may 
be expressed as ip = J2 X J2j i ) ( x -,j)\ x ) ® \j)- 

For a QW launched with initial state V'Oj the probability P(x, t) of finding the walker at 
the position x £ Zjv at time t is given by the standard formula P(x, t) = Tr (\x)(x\p(t)), 
where the time-dependent density operator p(t) is defined by: 

/J (t) = ^J = [/ i |Vo)(^o|^ t *. (1) 

During the evolution of the quantum walk, the history of decoherence-inducing events, 
if any, including any acts of measurement, may be modeled by the probabilistic option 
of applying to the coin degree of freedom, at each time step of the walk, a unital family 
of operators {A n }o< n < u , jointly satisfying the condition: 

E & n A n = I. (2) 

0<n<i/ 

Accordingly, when adjusted for decoherence, the density operator of the system acts 
on the probability density function p via the formula: 

P (t + i)= E uA nP (t)Altf. (3) 

0<n<v 

For simplicity, and without loss of generality, we may assume, from this point on- 
ward, that every quantum walk under consideration is launched from position |0) in coin 
state ipo- Upon applying Fourier transformations to all elements of the QW system, the 
formulation of the density operator assumes the following form: 

^) = ^EEiwi® £ ^°>^i' ( 4 ) 



k k' 



where |fc) = -4= J2 X Gxp(i2nxk/N)\x) and where the so-called "decoherence super-operator" 
Ckk' is defined by the formula: 

£w|Vo)^o| = E^WinlVoX^oli^AO (5) 

n 

A n {l,k,k') ^ 12 (l,fe,fc') 

A 2l {l,k,k') A 22 (l,k,k') 
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in which 



U c (k) 







2 77 / /,■ 

e at 



Cfc 



(6) 



Note that all essential features of the quantum walk, including all peculiarities of 
quantum behavior connoted by the words "decoherence" and "entanglement" arc fully 
encoded in the anatomy of the super-operator Ckk' ■ This is what enables us to study the 
level of quantum entanglement in quantum walks subject to decoherence. 



After t iterations, let the elements of the 2x2 matrix £ kk , |V , o)(V'o| be denoted by 



4fe'l^o)(^o| 



A n (t,k,k') A 12 (t,k,k') 
A 2 i(t,k,k') A 22 (t,k,k') 



(7) 



In terms of the standard basis for the 2iV-dimensional Hilbert space T-L = T-Ln <8> %, 
the density operator p(t) is given by: 

"(*)= E E p *vd*)\*)(v\®\3)(l\, (8) 

l<x,y<N \<j,l<2 

where 



p*vn(t) = ^J2J2( x \ k W\y)Mt,k,k'). 



(9) 



k fc' 



In terms of the above preliminaries, the probability P(x, t) of hnding the walker at 
position x at time t is given by: 

P(x,t) - Ti(\x){x\p(t)) 

iEE^i^^'i^^^i^)^!) 



N 



k fc' 



i ^ X>l*><*» (^ll(*. fc > ^') + 4m(*, fc, fc') 
fc fc' 

fa;xll(i) + ^xx22(i)- 



(10) 



To avoid unpleasant complications and to permit us more easily to illustrate our ap- 
proach to the analysis of a QW on the TV-cycle subject to decoherence-inducing influences, 
we concede, in this paper, to confine our attention to a specific model. Let /3 £ (0, : |) 
and fc <E {0, 1, ... , N — 1}. To serve as the coin operator of the system, we choose 



Ucifi) 

whose Fourier dual is given by 



cos j3 sin /3 
sin /? — cos j3 



e * » cos/3 e *» sin/3 
e^^sin/? — e^T^cos/? 



(11) 



(12) 



Note that when f) = j , the operator given by Eq. ([12")) is none other than the Hadamard 
coin operator. 

By the same token, to serve as the unital family {A„}o< n <i/ of decoherence-inducing 



Mathematical Structures in Computer Science 5 

operators on the coin degree of freedom, as in Eq.Q, we specialize to the following choice 
of three (v = 2) operators: 



A = V 1 -P°o, M = — (a + a z ), A 2 = —{o-q - <r x ), 

where < p < 1 and do and o~ z are the Pauli matrices. The level of decoherence induced 
by these operators is reflected by the value of p, called the decoherence rate. Specifically, 
the QW evolves as if the state of the coin is measured at each time step with probability 
p. Thus, when p = 0, the QW evolves as a purely coherent quantum process. At the other 
extreme, when p = 1, the QW behaves exactly like a classical random walk. 

Now let L(C 2 ) denote a the Hilbert space of all 2 x 2 complex matrices with inner 
product given by 

(M 1) M 2 ) = tr(M 1 t M 2 ). (13) 

Lemma 1. Let S be a superoperator on the Hilbert space L(C 2 ), defined by 

2 2 

S = J2 u ^ ■ A i u ? ■ B h> Y, UxAnBAlUz, 

ra=0 n=0 

where C/i, U 2 are 2 x 2 unitary matrices and B G L(C 2 ). Then (SB, SB) < (B,B). In 
particular, (SB, SB) = (B,B) for all B € L(C 2 ) if and only if the decoherence rate 
p = 0. 

Proof. See (|Liu and Petulante 2010| . 

An immediate corollary of this lemma, essential to our analysis, is the fact that |A| < 1 
for every eigenvalue A of iS. To justify this assertion, suppose that B\ is an eigenvector of 
S belonging to A. Then (SB\,SB\) = (XB\,XB\) — \X\ 2 (B\,B\). But since, according 
to the lemma, (SB, SB) <(B,B),we see that |A| < 1. 

Next, let us proceed to cast our reasoning in terms of the super-operator Ck,k'i which 
maps L(C 2 ) to L(C 2 ). If we choose as a basis for L(C 2 ) the Pauli matrices ao, a x , a y 
and o~ z , then, in terms of this basis, the 4x4 matrix representation of Ck,k' is given by: 



-k,k' 



c 


iqs~ sin 2/3 





is cos 2/3 





— qc + cos 2/3 


qs + 


c+ sin 2,5 





— qs + cos2/3 


—qc + 


s+ sin 2/3 


is~ 


qc~ sin 2/3 





cr cos 2/3 


ion 


we have set q = 1 


— p and 


,+ _ 


= cos \ >, S+ 


. 2t 

= sm — 


-(fc'+fe) 

N 



(14) 



_ 27r(fc'-fc) _ . 2ir(fe'-fe) 

C =COS ^ ' , s =sm ^ ; . 



After a somewhat tedious, but not very difficult calculation, we arrive at the the 
following explicit formula for the characteristic polynomial /(A) of Ck,k''- 
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/(A) = dct (AI 4 - C k , k .) 
= A 4 

+ (l + cos2£)(< ? c + -c-) A 3 
+ ((1 + q 2 ) cos 2/3 - 2qc + c-(l + cos 2/3)) A 2 
+<?(l + cos2,S) (c + -qcr) A 

V- (15) 

The following proposition summarizes some basic attributes of the eigenvalues of Ck,k'- 



Proposition 1. Suppose < p < 1. Let a typical eigenvalue of £fc,fc' be denoted by A. 
Then: 

(i) I|A||<1; 

(ii) If ||A|| = 1 then A = ±1; 
(iii) A = 1 when and only when k = k' , in which case the algebraic multiplicity of A = 1 

is 1 
(iv) A = — 1 when and only when |fc' — k\ = -y, in which case the algebraic multiplicity 

of A = -1 is 1. 

Proof. See Appendix A. 

Proposition 1 enables us to specify the long-term behavior of the matrix components 
of the total density operator given by Eq. © . 



Proposition 2. For the matrix defined in Eq. ([7]), the following assertions hold: 

(i) Suppose k = k' . If j = I then lim^oo Aji(t, fc, k') = =. If j 7^ I then lim^oo Aji(t, k, k') 

0. 
(ii) Suppose \k - k'\ = f . If j = I then ]im. t ^ 00 {-l) t A j i(t,k,k') = §. If 3 ^ I then 
lim t _ >00 Aji (£, fc, k') = 0. 

Y and 7^ 0. Then, for all combinations of 3, 1, we have limt^oo Aji 



(iii) Suppose \k— k'\ 7^ y and 7^ 0. Then, for all combinations of j, I, we have lim i _ i . 00 Aji (i, k, k') 



0. 
Proof. See Appendix B. 

In view of Proposition 2, the behavior, ast^ 00, of the operator in Eq. (jHJ can be 
specified as follows. 

Theorem 1. For a quantum walk on the iV-cycle, let the total density operator p(t) be 
defined as in Eq. ((8|). 

(i) Suppose N is odd. If x — y and j — I, then lim.t_j.oo Pxyji(t) = jft- If x 7^ y or j 7^ I, 
then limt-j.oo P xy ji(t) = 0. 
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(ii) Suppose N is even. If x = y and j — I and if t — x is even, then lim t _ i . 00 P X yji (t) = ■h. 

Otherwise, lim^oo P xy ji{t) = 0. 
Thus, in every case, as t — > oo, the operator p(t), viewed as a 2N x 2N matrix, converges 
to a diagonal matrix. If N is odd, the diagonal elements all converge to ^ . If N is even, 
then the diagonal elements converge, in an alternating pattern, to jr or 0. 

Proof. See Appendix C. 

As remarked by Zureck in (Zurck 2003), decoherence on a quantum system is mani- 
fested through its density matrix by the vanishing of the off-diagonal elements. In the 
context of a coin-driven quantum walk on the iV-cycle, this is precisely what Theorem 1 
asserts. Note that the off-diagonal elements of the density matrix are precisely the ele- 
ments which represent the quantum correlations (a.k.a. entanglement) between the coin 
subsystem and the position subsystem. Not surprisingly, at least for quantum walks, 
decoherence turns out to be practically synonymous with "disentanglement". A more 
precise elaboration of the relationship between decoherence and entanglement is deferred 
to the next section. 

For the position distribution P(x, t) (see Eq. (fTU| ). Theorem 1 has an immediate corol- 
lary, which echoes a similar result given in (jLiu an d Pctulantc 2010). 

Corollary 1 . Suppose a quantum walk on the iV-cycle is launched from the origin with 
initial coin state |V>o) and with decoherence rate p > 0. If N is odd, then P(x,t) con- 
verges to jf on all nodes of the cycle. If N is even, then P(x, t) converges to ■£ on the 
supporting nodes of the cycle and to on the non-supporting nodes of the cycle. 



3. Entanglement versus decoherence 

Since a purely coherent QW is a reversible process, the von Neumann entropy of the 
total density operator is invariant relative to time. If the QW is launched in a pure 
state, then it will continue to evolve in a pure state and the entropy of the reduced 
density operator on the coin subsystem can serve as a measure of its degree of entan- 
glement relative to the subsystem of the walker (jCarneiro et al. 20051 lAbal et al. 20061 
|Venegas-Andraca and Bose 2009[ TAnnabestani et al. 2010a|t . 

Unlike the purely coherent case, a quantum walk, when subject to decoherence, evolves 
in a mixture of states. If the "decoherence rate" is non-zero, then the von Neumann 
entropy of the total density operator is invariant no longer. Based on a very simple 
argument, it can be shown that the entropy of the total density operator is a strictly 
increasing function of time. The following two facts constitute a basis for the argument: 
(a) projective measurements increase entropy and (b) the entropy is a concave function 
of its inputs. Thus, to measure the level of quantum entanglement in a QW subject to 
non-zero decoherence, the von Neumann entropy must be considered separately for each 
of the subsystems as well for the total system. 

The von Neumann entropy of a quantum system A, denoted S(A), is a measure of the 
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uncertainty implied by the multitude of potential outcomes as reflected by its density 
matrix p(A). By definition, S(A) = S(p(A)) = -Tr(plnp). 

For a composite system with two components A and B, the joint entropy of their 
conjunction, denoted by S(A, B), is defined by the formula S(A, B) — — Tr(p AB lnp AB ), 
where p is the density matrix of the composite quantum system AB. 

A good measure of the level of quantum entanglement between the two components A 
and B is the so-called mutual information S(A : B), defined by the formula S(A : B) = 
S(A) + S(B)-S(A,B). 

We continue to confine our attention to the case of a quantum walk on the iV-cycle 
subject to decoherence on the coin degree of freedom. For the remainder of this article, 
our main objective is to show that whenever the decohrence rate p > 0, the mutual 
information between the subsystem of the coin and the subsystem of the walker eventually 
must diminish to 0. 

The following Lemma, due to Watrous (Watrous 2008), is essential to our reasoning. 

Lemma 2. Let X denote a complex Euclidean space and let Pos(A") denote the set 
of positive semidefmite operators defined on X with norm ||/o||tr = Tr I yp^p)- Then, 
with respect to this norm, the von Neumann entropy S(p) is continuous at every point 
pePos(X). 

Proof. See (| Watrous 2008^ . 

The operator norm || ■ || tr , which appears in the statement of Lemma 2, is known as 
the Schatten 1-norm or sometimes simply as the "trace norm" . 

Theorem 2. Suppose the QW is launched on the TV-cycle with initial coin state l^o) 
and with decoherence rate p > 0. let p(t) denote the time-dependent density operator 
of the overall system. If N is odd, then lim t _ i . 00 S(p(t)) = 1 + log 2 N. If N is even, then 
lim^oo S(p(t)) =\og 2 N. 

Proof. See Appendix D. 

Recall a basic fact from information theory: for any operator p defined on a Hilbert 
space of dimension d, the maximum value of the entropy S(p) is log 2 d. This implies that 
the limiting entropy values given by Theorem 2 actually are maximal, both for even and 
for odd values of the cycle length N. To see this, note that the overall Hilbert space 
% x WjV) over which the quantum walk evolves, is of dimension 2N. This easily explains 
why, when N is odd, the entropy is maximal. However, when TV is even, exactly half 
of the nodes of the cycle (every second one) are necessary and sufficient to determine 
completely the position distribution of the walker's itinerary. So, when N is even, the 
quantum walk evolves over a Hilbert space which effectively is of dimension N. 

Finally, we consider separately the long-term trend of the entropies associated with 
the reduced density operators of the two constituent subsystems (coin and position) and 
the mutual information between them. 

For the subsystem associated with the coin, the time-dependent reduced density oper- 
ator pc(t) is given by pc(t) = tracew(/o(£)), where the subscript w signifies exclusion or 
"tracing out", relative to the overall system density operator p(t), of the walker's degrees 
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of freedom. Similarly, for the subsystem associated with the walker, the time-dependent 
reduced density operator pw(t) is given by pw(t) — tracec(p(£)), where the subscript c 
signifies exclusion or "tracing out", relative to the overall system density operator p(t), 
of the coin's degrees of freedom. 

The following theorem summarizes our main finding. 

Theorem 3. Suppose a quantum walk is launched on the iV-cycle with initial coin state 
|^o) and with decoherence rate p > 0. Let pc{t) and pw(t) denote, respectively, the time- 
dependent reduced density operators associated with the subsystems of the coin and the 
walker. Then the long-term trend of the mutual information between the coin subsystem 
and the walker subsystem is given by lim^oo S {pc{t) ■ Pw(t)) = 0. 

Proof. See Appendix E. 

Finally, it can be shown, without much difficulty, that as t — > oo, the entropy values 
of the reduced density operators pc{t) and pw(t) each converges to its maximum value. 
It suffices to employ some elementary calculus based on Theorem 2, together with the 
following basic facts: (1) S(p c (t)) < 1, (2) S(p w (t)) < In TV when N is odd and (3) 
S(py/(t)) < In y w ben iV is even. 

Corollary 2. Suppose a quantum walk is launched on the TV-cycle with initial coin state 
|?/>o) and with decoherence rate p > 0. Then 
(i) lim^oo S{p c {t)) = 1. 

(ii) If N is odd, then limt-^ S(p w (t)) = InN. 
(iii) If N is even, then lim^oo S(pw{t)) = In ^. 

4. Conclusion and further questions 

The model of decoherence used in this article is only one of several prevalent in the current 
literature. It would be interesting to investigate how quantum entanglement responds to 
other models of decoherence, and not just for quantum walks on the iV-cycle, but for 
quantum walks over other kinds of topological networks as well. We speculate that, in 
every case, decoherence serves to erase quantum entanglement between the subsystems of 
a given quantum system. This implies that, as remarked by Zureck in fZurek 2003J, the 
progressive disappearance of entanglement should be accompanied by the correspond- 
ing disappearance of the off-diagonal elements of the density operators. Ultimately, the 
density operators should become indistinguishable from diagonal matrices. 

The natural question which arises in this context is: how long does it take for a quantum 
system, subject to decoherence, to reach its stationary state, devoid of any entanglement? 
To address this question, we propose a definition of decoherence time analogous to the 
measure of mixing time used in literature. Let p^ denote the limiting (stationary) density 
operator of the quantum system and let || • ||t r denote the trace norm as defined above. 
Then, for every e > 0, we define 

J D(e) = min{r|Vi>r: \\p(t) - Poo \\ tI < e} . (16) 

To estimate D{t), it would suffice to have good control over the eigenvalues and eigen- 
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vectors of the superoperator C l kk , . Unfortunately, even for very simple systems, the task 
of specifying the eigenvalues and eigenvectors of the superoperator L\ k , can be somewhat 
of a challenge. But the rewards of success would be worth the effort. Among other things, 
it would permit us to compare mixing time and decoherence time. The relationship could 
prove quite interesting. 
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Appendix A. Proof of Proposition 1 

Proof of (i). Ck.k' is a special case of the superoperator iS in Lemma 1, according to 
which, the moduli of all eigenvalues of Ck.k' are less than or equal to unity. 

Proof of (ii). Suppose e 10 is a non-real eigenvalue of Ck,k', where 9 is a real number. 
Then the conjugate e~ t6 also must be an eigenvalue and e~ t6 ^ e l9 . Hence /(A) = 
(A - e i6 ){\ - e~ ie )[\ 2 + aX + (1 - p) 2 ] for some a e C. For brevity of notation, let 
q = 1 — p. Comparing corresponding coefficients of both sides of Eq. (fi~5|). we obtain the 
following system of equations: 

a-2cos6» = (l + cos2/3)(< ? c+ - c") 
1 + q 2 - 2a cos 9 = -2qc + cr + (l + q 2 - 2qc + c~) cos 2^ 

a-2<7 2 cos6> = (1 + cos 2/3) (qc + - q 2 c~) . (17) 

After some elementary algebraic manipulations, we infer that 

2 / o^ 27r(fc' + fc) 2-K(k'-k) 
1 + q 2 = -q{\ + cos 2/3) cos — K — '- cos — y — , 

which is impossible since the modulus of the LHS is strictly greater than the modulus 
of the RHS. More precisely, note that the modulus of the RHS is strictly less than 2q 
which, in turn is less than 1 + q 2 . This contradiction implies that any unit eigenvalue of 
Ck,k' must be real. 

Proof of (iii). A = l is an eigenvalue of Ck,k' iff/(l) = (1 + cos 2/3) (1 — cos ^ N ~ )[1 + 
2(1 —p) cos * N — - + (1 — p) 2 } = 0, iff 1 — cos -^jj — " = 0j which implies k' = k. Moreover, 
since /'(l) = (1 — cos 2/3) [1 — (1 — p) 2 ] ^ 0, the algebraic multiplicity of A = 1 is 1. 

Proof of (iv). A= —1 is an eigenvalue of £fe : fe' iff/(— 1) = (l+cos2/3)(l+cos ^ N ~ )[1- 
2(1 - p) cos 2 " ( ^ +fc) + (1 - p) 2 ] = 0, iff 1 + cos 22^) = 0, which implies \k' - k\ = f . 
In this case, since /'(— 1) = (1 — cos2/3)[(l — p) 2 — 1] ^ 0, the algebraic multiplicity of 
A = — 1 is 1. 



Appendix B. Proof of Proposition 2 

Proof. Viewed as a 2 x 2 matrix, C\ k , \ipo)(tl>o\ is a linear combination of the Pauli 
matrices <to, <7 x , <r y and a z with corresponding weights given by 



W(t,k,k') = Cl 



k,k> 



a i 

Oil 

04 
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where the column vector \a.\, 0.2,0*3, a.i\ T represents |"0o) (V^o | with a.\ = i for all choices 
of |Vo). 

Proof of (i) . For brevity of notation, let 



q = 1 — p, c = cos ■ 



47rfc 



4nk 
~N~' 



When k — k', we have: 



where 



Therefore 



Qo = 



£k,k — 


1 
_ 



Qo _ 


? 


qc cos 2/3 

—qs cos 2/3 
g sin 2/3 


qs 

-qc 




c sin 2/3 

s sin 2/3 

cos 2,5 


r f - 

*-fc,fc — 


' 1 







? 



(18) 



(19) 



(20) 



in terms of which, the four components of £i u\tpo){' l Po\ can be expressed as follows: 
An = |+ei) ^22 = 5 — £i and A12 = A21 = €2, where ei and £2 are linear of combinations 
of the elements of the matrix Q l . By Proposition 2, the eigenvalues of Qo all possess 
moduli strictly less than 1. Therefore Qq — > as t — > 00, from which follows assertion (i) 
of the proposition. 

Proof of (ii). Similar to the reasoning above, when \k 



fc'| = y, we obtain: 



C l - 

J -k,k' — 



i-iy 








(21) 



where Qi is the 3x3 matrix adjoint to the leading entry c~ of the 4x4 matrix £k,k' 
in Eq. (fl"4"|) . The conclusion follows by a line of reasoning analogous to that used in the 
proof case (i). 

Proof of (iii). If |fc — k'\ ^ ^ and ^= 0, then, by Proposition 1, the modulus of every 
eigenvalue of £k,k' is strictly less than 1. Thus L\ k , — > as t — > 00, which suffices to 
conclude assertion (iii). 



Appendix C. Proof of Theorem 1 

Proof of (i). Suppose N is odd. For clarity of presentation, we proceed by considering 
three sub-cases: (a) x = y and j = I; (b) x ^ y and j = I; (c) j ^ L 

Proof of sub-case (i)(a). Suppose x = y and j — I. By Eq. ([9]), we have 



N Z P T 



i(t) 



k=k< 
JV-1 



: (fc-fc') 

t? Ajj(t,k,k') + 



E l-rrix {k-k' ) 
e w A,j(i,M) 



^%(t,fc,fc) 



fc=0 



e » A j:j (t,k,k'). 

k^k' 



(22) 
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By Proposition 2, If k = k' then the term Ajj(t, k, k) in (f22j) converges to \. If k 7^ fc', 
then the term Aj (t, k, k') in (f22j) converges to zero. So lim^oo P xy ji(t) = ^ . 
Proof of sub-case (i)(b). Suppose x ^ y and j = I. By Eq. (|9]), 

AT-1 , 

N 2 P xyjj (t) = J2 ^^^Anit^'V+Y, e 2 " ( ^ Vk 'A J3 (t 7 k,k'). (23) 

As in the prior case, referring to Proposition 2, we conclude that lim^oo P X yji(t) = 0. 

Proof of sub-case (i)(c). suppose j 7^ I. In this case, by Proposition 2, Aji(t, k, k 1 ) — > 
as t — > co. Since P xy ji(t) is a linear combination of Aji(t, k, k') for k, k' — 0, 1, 2, ..., JV— 1 
and j, Z = 1, 2, it follows that P xy ji(t) = -^ £] fc Sfc'( x l^)(k'l2/)Aji(A fc > fc') -> as i ^ 00. 

Proof of (ii). Suppose N is even. As above, we proceed by considering the same three 
sub-cases: (a) x = y and j — I; (b) x ^ y and j = I', (c) j ^ L 

Proof of sub-case (ii) (a) . Suppose x = y and j = 1. By Eq. ([9]), 

x > 27rix(k-k') , ^ > 2^ix(k — k') , 

N 2 P xxn (t) = £>— *» A jj (t,k,k')+ ]T e " A 33 (t,k,k') 

k=k> |fc_fc'| = f 

El-wix ik-k' ) 
e w A 3J (t,k,k') 

AT-1 

= ^%(t,fc,fc)+ E (-lf%(i,fc,fc') 

fc=0 |fe-fc'| = f 

e w %(t,fc,fc'). (24) 

By Proposition 2, the first and third summands of Eq. (|24|) converge respectively to 
Y and 0. Meanwhile, the second summand converges to (— l) x ~*y. It follows that 
lim^oo P X xjj{t) — jh when a; — t is even and lim^oo P XX jj(t) = when x — t is odd. 
Proof of sub-case(ii)(b). Suppose x 7^ y and j = /. By Eq. ©, 

N 2 P xyn {t) = ^e 2 ^T H, %(i,fc,fc')+ ]T e I ^ ii e" 2 ^ ±: %(t,fc,fc') 
fc=fc' |fe-fe'|=f 

e —%—e » Ajj(t,k,k') 

AT-1 , 

E2TTifc(3: — tj) x ^ lirik (x — y) ,, 

e— ^^%(t,fc ) fe)+ J] (-ire^^^%(t,fc,fe') 

fc=0 |fe-fe'| = f 

+ E e a " ( "5T ,fc ' ) ^ i (t,*,A')- (25) 

By Proposition 2, Ajj(t, k, k) converges to ^ for all fc, which implies that the first sum- 
mand of Eq. (|25p converges to as t — » 00. Meanwhile, in the third summand, every one 
of the terms Ajj(t, fc, fc') converges 0. It remains only to evaluate the second summand, 
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which can be reconstituted as 

H_p ^ e i-^ ( _ 1)M .. M)fc/) _ (26) 

ife-fc'i=f 

By Proposition 2, lim t ^ 00 (— l)*A,j(£, fc, fc') = i. Thus, the third summand also converges 
to as t — > oo. In summary, we conclude that limt_>. 00 P X yjj(t) = 0. 

Proof of case (ii)(c). Suppose j j^ I. The assertion is proved in the same way exactly 
as in the corresponding case when TV is odd. See sub-case (i)(c) above. 

Appendix D. Proof of Theorem 2 

Proof. Let p(oo) denote the diagonal matrix diag(2^-, t^-, ■••, 2]\f)- If N is odd, Theorem 
1 implies that lim t _ ) . co ||p(i) — /?(oo)||tr = 0. By Lemma 2, it follows that lim^oo S(p(t)) = 
S(p(oo)) = 1 + IniV. When TV is even, the assertion lim^oo S(p(t)) — In AT is similarly 
justified. 

Appendix E. Proof of Theorem 3 

Proof. When combined with Theorem 2, the inequality S (pc(t) , pw (t)) < S(pc(t)) + 
S(pw(t)) < 1+lnAT implies that lim t _ KX) S(pc(t),/»w(*)) = 1+lniV. Thus, Um t ^. co [S'( / 9c(i))- 
S(pw(t))} = 1 + InJV. Since S(p c (t) : p w (t)) = S(p c (t)) + S(p w (t)) - S(p c (t),Pw(t)), 
therefore lim^oo S(H N ■ "He) = lim^oo S(p c (t) : Pw(t)) = lim t _ > . 0O [5'(pc(*))+S , (/-w(*))]- 
lim^oo S(p c (t), pw(t)) = 0. 



